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Abstract 

A representation of the continuum fermionic propagator as a sum 
of directed random walks on a lattice is presented. Also a random 
walk representation for the lattice fermionic propagators is developed 
for the case of the naive, the Wilson, and the Kogut-Susskind fermions. 

For the naive fermions the phenomenon of fermion doubling appears 
as having 2 ^ distinct spin factors being associated with a single path 
in D-dimensions. In the case of the Wilson and the Kogut-Susskind 
fermions, in the naive continuum limit, the path integral representa¬ 
tion coincides with the path integral representation for the continuum 
fermionic propagator. Using this representation the Green’s functions 
of lattice QCD involving quark operators are written as a sum over 
the paths of valence quark, the gauge fields and the sea quarks be¬ 
ing integrated out. Possible advantages of such a representation are 
illustrated by showing how one can use numerical simulations to ob¬ 
tain a heuristic insight into the relationship between QCD and the 
constituent quark model. 



1 Introduction 


It seems increasingly likely that in the near future the lattice simulations 
will provide an accurate numerical solution of Quantum Chromodynamics 
(QCD). Complimenting these numerical investigations are many attempts at 
obtaining a heuristic understanding of the dynamics of QCD by trying to 
isolate those fluctuations in the gauge fields that may be responsible for the 
conhnement and chiral symmetry breaking ||I[. A similar interplay between 
numerical simulations and heuristic understanding is lacking for the problems 
involving the quark degrees of freedom. Thus there are many striking features 
of strong interactions, an incomplete list would include the apparent successes 
of the non-relativistic quark model and the Okubo-Zweig-Iizuka rule [Q], for 
which the numerical simulations provide little or no heuristic guidance. 

The reasons for this are familiar, the fermionic degrees of freedom be¬ 
ing elements of Grassmann algebra cannot be simulated numerically. These 
degrees of freedom have to be integrated out from the functional integral 
before it can be numerically evaluated. The result of the integration is ex¬ 
pressed in terms of the fermionic determinant and the fermionic propaga¬ 
tor i i B 1, 0- But the numerical calculation of the determinant and the 
inverse of a fermionic matrix seems to offer little scope for visualizing the 
degrees of freedom that gave rise to them in the hrst place. 

In this paper I suggest an alternate way of dealing with the fermionic prop¬ 
agator which is more amenable to heuristic investigations. The approach is 
based on the path integral representation, or the equivalent random walk rep¬ 
resent at ion[|, of the fermionic propagator [^, |^, TO, n]. One way of obtaining 


a random walk representation for the lattice fermions is to expand the prop¬ 
agator for the Wilson fermions in a Neumann series BBi- An alternate 
approach, which is more general, is based on Polyakov’s observation that 
the path integral for a fermion can be regarded as a natural generalization of 
the path integral for a scalar particle by including an appropriate spin factor. 
In this paper I will take the latter approach. 

With in such an approach it is natural to ask the question that how does 
the phenomenon of fermion doubling on a lattice |T3, |T^ manifests itself in 


the language of path integrals. To answer this, in Sec. 2, hrst I represent the 


^In this paper I will interchangeably use the words path-integral and random walk 
representation. 
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continuum fermionic propagator as a sum over directed random walks on a 
lattice. Then in Sec. 3 and 4, I consider the path integral representation for 
the lattice fermionic propagators, this includes the so called naive fermions, 
the Wilson fermions, and the Kogut-Susskind fermions. In the case of the 
naive fermions the expected phenomenon of fermion doubling appears as the 
existence of 2^ distinct spin-factors for a single path in a D-dimensional hy- 
percubic lattice. Continuing the analysis in Sec. 4, it is shown that for the 
case of the Wilson and the Kogut-Susskind fermions, in the naive continuum 
limit, only one of the 2^ spin-factor survives and the path integral repre¬ 
sentation of the propagators coincides with that of the continuum fermionic 
propagator developed in Sec. 2. 

Having developed a path integral representation for a free fermionic prop¬ 
agator it is easy to extend it to the case when a gauge field is dehned on a 
lattice 0,0. Using the resulting representation of the fermionic propaga¬ 
tor, in Sec. 5,1 write the Green’s functions involving quark fields as a sum over 
the paths of the valence quarks. This way of writing the Green’s functions 
in which the integration over the gauge held is carried out before summing 
over the quark paths is a close transcription, which allows for numerical ex¬ 
plorations, of the heuristic way in which Wilson [|^ motivated his criterion 
for the quark conhnement. Perhaps the advantage of this approach is that 
it allows us to probe the theory using a language, the paths of the valence 
quarks, which is easy to visualize. This is illustrated in Sec. 5 by indicating 
how the numerical simulation can provide a insight into the relationship be¬ 
tween QGD and the constituent quark model. I state my conclusions in the 
last section. 


2 Representing Continuum Fermionic Prop¬ 
agator as Random Walks on a Lattice 


The fermionic propagator can be represented as a path integral either us¬ 
ing an appropriate spin-factor g, or equivalently as a sum over directed 
random walks [§, 0 0. In this section I will represent the continuum Eu¬ 


clidean Fermionic propagator as a sum over directed random walks on a D 
dimensional hypercubic lattice. Following 0 


I start with the Fermionic 
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propagator in the momentum space 


A(p, m) 


1 


( 1 ) 


where 7 ^ are the Euclidean Dirac matrices satisfying the anticommutation 
relation 

{ 7 /^, 7 ^^} =-2 ( 2 ) 

The propagator is now written as a Laplace transform 


1 /*oo jYf J\ _ 

A{p,m) = —J^ exp{— —}G{L,p)dL, (3) 

G(i.p) = exp{^^^}, (4) 

where D is the dimension of the Euclidean space which, unless stated other¬ 
wise, will be taken to be four. The motivation for expressing the propagator 
as a Laplace transform is that the Fourier transform of Eq. can be then 
interpreted as the probability amplitude for the particle to reach the point 
X by traveling along paths of length L. With this in mind, the Laplace 
transform is expressed as 


G{L,p) 


lim 

N—^oo 


DN 


N 


(5) 


The individual terms appearing in Eq. can be written as 


l-pL 

ND 


I 2D ^ 

^£(1 - - irV.p—), 


( 6 ) 


where fp are the directions available to a particle on a D-dimensional hyper- 
cubic Euclidean lattice, and for D = A can be taken to be: 

= (1,0,0,0), = (-1,0,0,0), 

= ( 0 , 1 , 0 , 0 ), = ( 0 ,- 1 , 0 , 0 ), 

= (0,0,1,0), = (0,0,-1,0), ^ ^ 

= (0,0,0,1), h® = (0,0,0,-1). 

The identity, Eq. (^, can be verihed by using the above choice of hh 
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This should be contrasted with the following representation 

f . L 


1 - 


ND 


/ I ^ 

dh{l — i'y.h){l — ih.p—) 


( 8 ) 


where the integration is over all the directions in D-dimensional Euclidean 
space and results in a representation of the continuum fermionic propagator 
as a sum over paths in the continuum. 

Substituting Eq. @ into Eq. (|]) leads to 


N 


2D 


G{L,p) = n 


k=l 


i=l 


N' 


(9) 


which in the limit iV —> cx), for a fixed L, can be writen as 


N 2D 


2D 


G{L,p) = (—) ^ 

i^ = \ i\=\ 


L ^ 

X exp(-z—p. 

k=l 


( 10 ) 


Taking the Fourier transform of Eq. (IT) leads to 

1 V 2D 2D 

G{L,x-y) = E ■■■ 


= l 2l = l 

L ^ 


x6^(x — y - n 

k=l 




( 11 ) 


This, as anticipated, can be interpreted as a sum over paths of length L 
starting at the point y and ending at the point x. The paths in Eq. ( 0 ) are 
defined on a D dimensional Euclidean lattice, with each step in the direction 
h* there is a spin-factor (1 — ij.fp) associated with it. The step size defines 
the lattice constant 

In view of Eq. (|T^ the fixed length propagator can be written as 
G{L,x-y) = G{N-,a-,x -y), 


N 2D 


2D 


G{N-,a-,x-y) = bm ( —) ^ ^ (1 - • (1 - zq.h*!) 

ijv—1 ii=l 

u [X — (/, u 2_ 

X 


-0 V2D/ 

6^{x - y,aY.k=i b*'') 


a 


D 


(13) 


4 







In Eq. ( p!3[) the Dirac-h function has been replaced with the Kronecker-J 
function. 

Using Eq. and replacing the integration over L by a summation over 
N leads to the following expression for the continuum fermionic propagator 
in the coordinate space 


A(x -y,m) = lim — ^ exp( 


N=1 


maN 

D 


)G{N; a;x-y). 


(14) 


Though the sum in the above equation starts from N = 1 but in the limit 
a —>■ 0, which is equivalent to the \\ x — y ||>> a, only terms with very large 
value of N contributes because of the presence the 5-function in Eq. ([TB|). 
Above representation of the continuum fermionic propagator as a sum over 
directed random walks on a lattice can be extended to the case when a gauge 
held is dehned on a lattice. Following 0. see also [|T^, the continuum 


fermionic propagator in the presence of a background gauge held, U, which 
is dehned on a lattice can be written as: 


A(x - y] m; U) = hm ^ X! ®^P(' 


maN 

D 


)G{N;a;x-y,U), (15) 


where G{N; a;x — y; U) is given by 


G{N;a;x — y;U) = 


1 N 2D 

= ("-L^ 


2D 




H H (1 - • (1 - 

2i=l 


UfilJ „ (fp) 


a 


where U^k is a SU{N) matrix associated with the link of the lattice that 
corresponds to the kt^ step of the random walk. 

It is interesting to note that in the quenched approximation, where one 
neglects the contribution of the fermionic determinant, using above represen¬ 
tation of the fermion propagator one can give a non-perturbative dehnition of 
QCD that keeps the fermionic degrees of freedom in continuum while dehn- 
ing the gauge degrees of freedom on a lattice. This would be of signihcance if 
one could also dehne the fermionic determinant in a similar manner (for some 
attempts along this direction see f^), for this would give a non-perturbative 
dehnition of QCD that avoids the problem of fermion doubling. 
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3 Random Walk Representation of the Naive 
Lattice Fermionic Propagator 


In this and the next section a representation of the lattice fermionic propa¬ 
gators as a sum over paths on a lattice will be obtained. The naive lattice 
fermionic propagator 0.@.0i in the momentum space is given by 


An{p,m) 


1 

m + 


(17) 


where a is the lattice constant and the y^s are the Euclidean Dirac matrices 
satisfying Eq. (^. It is well known that in the continuum limit Eq. (0) 
describes the propagation of 2^ distinct fermions fl^, an example of the 
general phenomenon of fermion doubling |]^|. Writing the naive fermionic 


propagator as a sum over paths will allow us to understand this phenomenon 
in the language of path integrals. 

The development of the path integral representation for the naive lattice 
fermionic propagator starts, as before, with a formal Laplace transform of 
Eq. (|^ 


An{p,m) 

- ^ / exp{ }Gn{L,p)dL, 

JJ Jo D 

(18) 

Gn{L,p) 

-7^sin(p^a)L 

= “Pf Da > 

(19) 

and with the identity 



7 ^sin(p^a)L 

DNa 

1 2D ^ 

(20) 


where h* are the unit vectors, along both the positive and negative directions 
of a D-dimensional Euclidean hypercubic lattice and are dehned in Eq. (|^). 
Using Eq. (|T9|) and Eq. (^0]), as in section 2, we can write Gn{L,p) as 

T N 2D 2D 

Gn{L,p) = lim (—) ^ 

N^oo ZU ■ \ 

lN = l 2l=l 

L ^ 

xexp(-i—sin(p^a).^h;'=)}. (21) 
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The Fourier transform of Eq. (pT]), defined by 

■f d^p 


GniL,x- y) = 


exp{ip.{x - y))GniL,p) 


'dr (2vr)^ 

with X and y belonging to the lattice, can be written as 

N 2D 2D 

Gn{L,x-y) = ^ ■ x;(l - "'(I - 


( 22 ) 


ijY=l ii=l 

x^lx-u, 


(23) 


where I[x — y, {h*}] is given by 
'f d^p 


I[x-y,{n^}] = 


L 


^ (27r) 


■^exp{ip.{x-y))exp{-i—sm{pf,a).J2f>''i^}- (24) 


N 

c 

k=l 


Next consider Gn{L, x — y) in the limit a —>■ 0. For this it will be convenient 
to, using the periodicity of the lattice momentum space, shift the range of 
integration from ^ ^ to ^ < |^. Also the range of for each p^ 

is divided into two regions, in the first region p^ ranges from ^ < P /2 < 
and in the second region it ranges from ^ < P/t < |y- In fhis manner the 
lattice momentum space divides into 2^ hypercubes centered around points 
for which p^, the pth component, is either zero or These centers will 
be denoted by the vector K^. Using this division of the momentum space 
Eq. can be written as 


20 

I[x-y]{fd}] = Y,h[x-y]{fd}], (25) 

h=l 

with Ih[x — p; {fd}] given by 

f d^p L ^ 

Ih[x - p; {h*}] = exp{-i— sin(p^a) ^ exp{ip.{x - p)). (26) 

In Eq. (^), as noted above, the pth component of p^ lies either in the range 
^ < P /2 < ^ or in the range ^ ^ P^ < |y- Dehning a new integration 
variable 

PkL = + ^/ 2 ! ( 27 ) 
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allows one to write Eq. 
J = 


as 


2a 


i2< 


D 


, L sm{Kl^a + k^a) ^ ^ 

exp(-^—-^ n^) exp{ik.{x - y)). 

(28) 


N 


k=l 


Taking the limit a —>■ 0 of the above expression leads to 


hni4[a;-2/];{h*}] = - y) - ^ cos(irJa)/c^. ^ h;"), (29) 


L 


N 


a^O 


N 


k=l 


where cos(K^a) is either +1 or —1, depending on whether is 0 or Thus 
Gn{L,x) in the limit a —0 has the following form 

-I N 2D 2D 

GniL,x) = (—) ^ ^(l-i7.h*^)---(l-z7.h*i) 


^2D 

2° 


i^ = \ 2l=l 


X y^5^{x - — cos{Kla)k^. ^ (30) 


h=l 


N 

E 

k=l 


Gn{L,x) can now be interpreted as a sum over 2^ distinct random walks of 
length L but for the factor of cos{K^). This factor can removed from the 
6 function by suitable redefinition of the spin-factors. To see this consider 
the case when the uth component of is then substitute the set of unit 
vectors {h*} by a new set of vectors obtained by changing the sign of the Uth 
component of the vectors {h®}. This transformation keeps Eq. (^OD invariant 
as it corresponds merely to relabeling of the vectors {h®}. Because of this 
the spin factor for a link 1 — i(- ■ ■ -1- 7 jyh(, + ■••)) where no summation over 
1 / is implied, changes to 1 — z(- • • — 7 i/h(, -[-■■■). This can be brought to its 
original form by doing a similarity transformation on the gamma matrices 



II 

1 

= 5^(5'®)-^ 

(31) 

result of these transformations 

Gn{L, x) can be written as 


Gn{L,x) 

2^ 

= J2^hiL,x), 


(32) 


h=l 




1 N 2D 

2D 


Gh{L,x) 

= (^) S- 

^(l-^7^h®^)•••(l-i7^.h®®) 



2JV = 1 

21 = 1 



X 

I 

h®'=), 

(33) 
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which can be interpreted as a snni over paths of length L with the step-size a 
where a is the lattice constant. Finally writing the formal Laplace transform, 
Eq. (P!B[), as a snm over lengths in lattice nnits leads to 


, a ^ . maN ^ ^ ^ 

hmA„(x,m) = 2 ^ exp(- —)}^Gh{N,x), 

_ T N 2D 2D 

G,{N,x) = 

iiv=i n=i 

aX> 

The above eqnation shows that in the continnnm limit the naive fermionic 
propagator represents propagation of 2^ distinct fermions. In the same limit 
the propagator of each of these fermions can be written as a snm over paths, 
or random walks, on the lattice. Different species differ by the spin-factor 
associated with them and by an overall phase factor. Since the phase factor 
can be absorbed in the normalization of the single particle wave fnnction and 
the spin factors differ only by a similarity transformation of the 7 -matrices, 
therefore we have in Eq. (|3^ a simnltaneous propagation of 2^ Dirac parti¬ 
cles. Farther analysis along the line of Ref. [Q revels that these donblers 
appear in pairs with opposite chirality. 


(34) 


4 Random Walk Representation for Wilson 
and Kogut-Susskind Fermions 

Having elncidated the phenomenon of fermion donbling in the langnage of 
path integrals for fermions, it is natnral to see how the Wilson and the Kogut- 
Susskind fermions avoid or mitigate this problem. I start with the Wilson 
fermions, the corresponding propagator in the momentum space is i, i 0 

^W{P) = 1 . . - ^ n - 1 - 

m + -sm(p^a) + - - cos(p^a)) 

The development of a path integral representation for the Wilson fermions 
is similar to the previous two cases, again one writes the propagator as a 
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formal Laplace transform 

Aiy(p) 


1 

D Jo 


dLexp{-^)GwiL,p), 
I JJ 

T D 


~ L ^ 

Gw = exp( —( 7 ^sin(p^a) + “ cos(p^a))). 

n=l 


(36) 

(37) 


Then using the following identity 
L ^ 

1 - + 11(1 - cos(p^a))) = 

M —-1 


1 2D 

1=1 

which can be verihed 
GwiL,p) as 


n.ft‘)(i - 

with fh 


(38) 

■ cos(ppa))}, 


'* as the unit vectors defined by Eq. , to write 


1 AT 2D 2D 

GwiL,p) = ji]^(—) ^ ■■■ (39) 

1 * 1=1 


^2D' ^ 

iiv=l * 1=1 


L ^ ■ L ^ 

X Sin(p^a). II ^ I] (1 - cos(p^a))). 

Taking the Fourier transform of Eq. (P9[ ) leads to 

1 TV 2D 2D 

Gw{L,x) = Ji^(—) H Il(l-Z7.h*^)---(1-Z7.h*i) 

1 21=1 


N—^oo 

xlw 


2D^ ,1 

2iV = l 

cp{ip.x) exp{—ij^ 


N 


(40) 


ex' 


sin(p^a). I^h;, 


Ak 


k=l 


Iw[x,{ff}] = 

L ^ 

- 75 -II(l-cos(p^a))). 

/*=! 

Eq. (IT) differs from Eq. (^) only by the presence of the Wilson term 
fore one can again write it as 

Iw[x] {fd}] = H Iwh[x] {h*}], (42) 

h=l 


(41) 

there- 
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where Iwh[x] {h*}] is given by 


Iwh[x-,{h"}] = exp(iA;.x) 


X exp —^ 


D 


(27r)^ 

. L sin(ir^a + k^a) ^ 

« ti " 




- ^(1 - cos(Ar^a + /c/,a)) 
M=i 


( 43 ) 


In the limit a —> 0 one obtains 


\\m lwh[x] [n"}] = e*-^'‘-^exp(-^ ^(1 - cos(A:Ja))) 


D 


a^O 




^ cos(AT^a) ^ h;'=). 

fc=l 


TV 


(44) 


From Eq. (|4^) one sees that in the limit a —> 0 only Iwh corresponding to 
= 0 will survive. Rest of the Iwh goes to zero, approximately as exp(—■^), 
and therefore their is no contribution from the “doublers” and one can write 


L ^ 

lim/iy[a:;{h*}] = 5^(x-— ^h*). (45) 

k=l 

Substituting Eq. (^5]) in Eq. (|40|) leads to the identihcation of L as the length 
of the path and N as the number of steps of the random walk. Finally the 
propagator in the continuum limit can be written as 


lim Avi/(x, m) 

a^O 


G\y{N', CL] x) 


n ^ rri n ^ 

— exp(- —)GwiN]a]x), 

^ N=1 ^ 

-i N 2D 2D 

(fT^) E ■ ■ • • • • (1 - 

ijv=l n=l 

^^5^(x-?/,aEf=ih*'=) 

X- ^ -. 


(46) 


(47) 


Thus, in the continuum limit one can represent the propagator of the 
Wilson fermions as a sum over paths on a lattice, further more the rep¬ 
resentation is identical to the random walk representation of the continuum 
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fermionic propagator. A path integral representation for the Wilson fermions 
in the form of a hopping-parameter expansion is well known and is similar 
to the path integral representation developed above. But conceptually the 
two differ, in obtaining Eq. (^) no expansion in hopping parameter is done 
and the path integral representation emerges only in the continuum limit. 
Conceptually the path integral representation for the Wilson fermions is a 
way of transcribing the spin factor, present in the continuum path integral 
representation of a fermion, on to a lattice . 


Next I consider the propagator for the Kogut-Susskind fermions. Since 
the development is almost identical to the previous cases only the hnal results 
will be stated. The propagator for the Kogut-Susskind fermions P, can be 
written as 


^Ks{p) = 


m + ( 7 ^ (g) l)i sin(p^ 6 ) + i E 2 =i (1 - cos(p^ 6))(75 0 ’ 


(48) 


where are 212 ] dimensional flavor matrices, and are given by 

t5 = 75, (49) 

and b = 2a, a being the lattice spacing. The propagator, Eq. (|^), describes 
the propagation of 2['2'1 degenerate flavors. Formally the Kogut-Susskind 
propagator differs from the Wilson propagator only by the presence of the 
flavor space matrices, and that the lattice spacing has been effectively dou¬ 
bled, so one can immediately write down the path integral representation for 
the Kogut-Susskind propagator as: 


\imAKs{x,m) 

a^O 


b 

D 


£ exp(-^^(^)G'xs(A^; b;x) 

N=1 ^ 


(50) 


where Gks is given by 


GKs{N;h;x) = 


1 N 2D 

(^) 


2D 


• yO ((1 — ■ (1 — zy.h*^)) ® 1 


ijv=l 


ii=l 


X 


b^ 


12 


(51) 







In Eq. one has the path integral representation for the propagator of 
2 [t 1 degenerate fermions. The paths are dehned on a blocked lattice with a 
step size of b = 2a. 

The path integral representation of the Kogut-Susskind propagator is, 
apart from the fact that it represents the propagation of 2 [' 2 l degenerate 
fermions, identical to the path integral representation of the continnum 
fermions. The fact that the path integral representation of both the Wilson 
and the Kogut-Susskind fermions, which are derived from different lattice 
actions and have different lattice symmetries, coincides with each other is 
not surprising as these representations appears only in the continuum limit. 
What is, perhaps, surprising is the fact that the continuum fermionic prop¬ 
agator can be represented as a sum over paths even when the paths are 
restricted to a lattice. 


5 Path Integral Representation of Green’s 
Function in QCD and the Quark Model 

As stated in the introduction one of the motivations for developing a path 
integral representation for the lattice fermionic propagators was to look for 
a formalism that can provide some heuristic insight into the role of quark 
degrees of freedom in QCD. To this end, in this section, I will write the 
Green’s functions of lattice QCD involving quark degrees of freedom as a 
sum over the paths of valence quarks thus making a connection with the 
heuristic language of quark model. 

Consider a meson propagator build from '^(a;)r'0(a;) as the interpolating 
held for a quark-antiquark meson with T determining the spin, havor, and 
the parity of the meson i0- Such a propagator can be written as 

{i>{x)Tij{x)^{y)Tij{y)) = Z-^ J DUeMSeffiU)) 

x(Tr[A(x,y,U)rA(y,x,U)r] 

-Tr[A(x, X, U)r]Tr[A(y, y, U)r]) (52) 

where the quark degrees of freedom have been integrated out giving rise to 
an effective gauge action and fermionic propagators. The partition function. 
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Z, and the effective gauge action, Seff{U), being deffned by 

2 = 1 DUexp(-Seff(U)). 
SeffiU) = +lnDetA-"((/). 


( 63 ) 

(54) 


Sg{U) represents the lattice regularized action for the gauge field U il- 
For the sake of definiteness, the fermionic determinant is assumed to be 
defined for the Wilson fermions. In what follows it will be convenient to use 
the following notation to represent the functional integration over the gauge 
fields 

{f{U))u = Z-^ j DUeM-Seff{U))f{U). (55) 

Also the path integral representation of the quark propagator, Eq. ([T5|), will 
be written in the following compact notation 


A{x,y,U) = ^exp(-/i/)<h(/a.j/)f/(^,,y), 

Ixv 


(56) 


where Ixy denotes a path on a lattice of length /, in lattice units, starting at 
y and ending at x, /u is a measure of the bare mass and is given by 


h = 


ma 


(57) 


while the spin factor ^{Ixy) and the gauge field factor U{lxy) are given ex¬ 


can be written as 


plicitly by Eq. 

In the above deffned notations Eq. 


{ilj{x)Vil){x)ij{y)V^lj{y)) = (Tr[A(x, y, U)rA(y, x, U)r])u 

-(Tr[A(x, X, U)r]Tr[A(y, y, U)r])u . (58) 


Consider the two terms appearing on the right hand side of Eq. (|58D sepa¬ 
rately and replace the fermionic propagators appearing in them by their path 
integral representation Eq. (|56|). The first term can be written as 


(Tr[A(x,y,U)rA(y,x,U)r])u = + lyx)) 

^xy i^yx 

xTr[<I)(hy)r<F(ly,)r] 

x(Tr(U(hy)U(lyJ))u. 


(59) 
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Similarly the second term can be written as 


(Tr[A(x,x,U)r]Tr[A(y,y,U)r])u. = + lyy) 


xTr[<l>(U)r]Tr[<h(l,y)r] 
x(TrU(U)TrU(lyy))u. (60) 

Another quantity of interest is the order parameter for the chiral symme¬ 
try, '0-0, which can be written as a sum over closed paths 


(^(0)^(0)) = (Tr[A(0,0,U)])u 

= ^exp(-/i/oo)Tr[$(loo)](TrU(loo))u- 


(61) 


^00 


In a similar manner all the Green’s functions involving quark helds can 
be written as a sum over one or more closed (valence) quark paths, each path 
being weighted by its length, its spin factor, and the expectation value of the 
path ordered product of the gauge helds along the quark paths. The effect of 
the dynamical or the sea quarks being included in the expectation value of 
the gauge held factor. The advantage of writing the quark Green’s function 
using the path integral representation of the quark propagator is that it 
separates the contribution coming from the motion of the valence quark, in a 
sense the kinematics, from the contribution coming from the huctuations in 
the gauge helds and the quark degrees of freedom. As a possible application 
of the random walk representation of the fermionic propagators, consider 
the relationship between QGD and the constituent quark model [^. In the 
chiral limit the constituent quark picture for pion is not useful, pion being 
the Nambu-Goldstone boson, while it is a good approximation for the rho 
meson. One would like to understand, starting from QGD, the emergence of 
the constituent quark picture for the rho and its absence for the pion. For 
this purpose consider the connected part of the meson propagator Eq. (|5^), 
a priori all possible paths contribute to this propagator but one expects 
that because of the conhnement the paths in which quark and anti-quark are 
widely separated are strongly suppressed. This motivates us to restrict the 


sum in Eq. (M) to a class of paths in which the quark and the anti-quark are 
separated at most by, say rmaxi where rmax is some measure of the meson size. 
One expects that the constituent quark model should emerge at a coarser 


level of description |21] and such a coarser description should appear as a 
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result of a renormalization group transformation. In the present context one 
could try and see if the constituent quark model emerges after the paths have 
been blocked, where blocking of paths is in analogy with the blocking of spins 
under the real space renormalization group transformation . One possible 


way of blocking the paths is by summing over the paths in the neighborhood 
of a given path belonging to the above dehned class. In particular consider 
a path Ixyx (where l^yx being a compact notation for the closed path Ixylyx), 
imagine that that this path is enclosed in a tube of radius Vb ~ ‘2rmax, and 
then sum over all the paths which lie inside this tube to create a blocked 
path. Denoting the paths lying inside the tube by I, one can conjecture that 
in the case of the rho meson the sum over these paths, 


^ exp(-;t<i)'iv|'i>(4)rp4(y )r,] {'ivu(4j>u, 


( 62 ) 


can be approximated by a single path with the following weight factors 


exp{-meff{lxy + /ya:)Tr[$(hy)rp<I)(lyx)rp] exp(-Vefr(lxy , lyx)) (63) 

where nie// is a measure of the constituent quark mass and exp(—is the 
weight factor coming from the conhning potential of the constituent quark 
model. The simplest possible guess for 14// would be aAxyx where a is the 
string tension and Axyx is the minimal area enclosed by the closed path Ixyx- 
On the other hand for pions one would conjecture that the blocked paths 
should be approximated by a single effective path with the following weight 
factors exp{—m.^lxy), where being the measure of the pion mass on the 
lattice which in turn should be related to the expectation value of the order 
parameter for the chiral symmetry breaking, ('44) • Such heuristic conjectures 
can be numerically tested, for the sum in Eq. (|6^) is well dehned on a lattice 
and is a small subset of all possible paths connecting point x and y. Perhaps 
the biggest obstacle in such a numerical investigation is the requirement 
of large lattice sizes, for the random walk representation appears only as 
one approaches the continuum limit. For hnite lattices one will have to 
understand how to separate the lattice artifacts from the continuum physics. 
A hrst step in such a direction would be to use the above formalism for 
exploring the dynamic of lattice QCD in two dimensions. 
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6 Conclusions 


Path integral or the random walk representation for the fermionic propagator 
was developed with two motivations in mind. One to try and understand the 
phenomena of fermion doubling on a lattice in the language of sum over 
paths, and other to look for a formalism for the fermionic degrees of freedom 
that might be more susceptible to heuristic investigations. 

In the first respect, it was shown that a naive lattice fermionic propagator 
can be written as a sum over paths and with each path there are associated 
not one but 2^ spin factors in D dimensions. As a result the naive propa¬ 
gator represents the propagation of 2^ fermions. This is in contrast to the 
case of the continuum fermionic propagator which can be represented as a 
sum over paths on a lattice with a unique spin factor associated with each 
path. In the case of the Wilson and the Kogut-Susskind fermions the path 
integral representation emerges in the continuum limit and coincides with 
the path integral representation of the continuum fermionic propagator. In 
this manner the path integral representation provides an additional insight 
into the relationship between various lattice fermionic propagators and their 
continuum counterpart, this may be of some use in exploring new ways of 
representing fermionic degrees of freedom on a lattice. 

As to the second motivation, it was shown that using path integral rep¬ 
resentation for the quark propagators, the Green’s functions of QCD can be 
written as a sum over the paths of the valence quarks. Such an representation 
of the Green’s functions allows for the possibility of delineating the important 
paths that contribute to a Green’s function. Also it allows for the possibility 
of implementing the ideas of real space renormalization group on variables 
which are easy to visualize and are close to the heuristic description of QGD, 
namely the paths of the valence quarks. This was illustrated by showing how 
one can use numerical simulations to understand the relationship between 
QGD and the constituent quark model. 

The path integral representation of the lattice fermionic propagator is 
unlikely to be useful for an accurate numerical solution of lattice QGD but 
it can be a useful tool for testing heuristic insights and for developing new 
intuitions. 
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